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ERROR ANALYSIS IN THE MEASUREMENT 
OF AVERAGE POWER WITH APPLICATION TO 
SWITCHING CONTROLLERS 

James E. Maisel 
Cleveland State University 

ABSTRACT 

Power measurement errors due to the bandwidth of a power 
meter and the sampling of the input voltage and current of a 
power meter were investigated assuming sinusoidal excitation 
and periodic signals generated by a model of a simple DC/ DC 
chopper system. Errors incurred in measuring power using a 
microcomputer with limited data storage were also considered. 

In this investigation, the behavior of the power measure- 
ment error due to the frequency responses of first order trans- 
fer functions between the input sinusoidal voltage, input sinu- 
soidal current and the signal multiplier was studied. Results 
indicate that this power measurement error can be minimized if 
the frequency responses of the first order transfer functions 
are identical. 

The power error analysis was extended to include the power 
measurement error for a model of a simple DC/DC chopper system 
with a DC power source and an ideal shunt motor acting as an 
electrical load for the chopper. The behavior of the power 
measurement error was determined as a function of the chopper's 
duty cycle and back emf of the shunt motor. Results indicate 
that the error is large when the duty cycle or back emf is small. 

Theoretical and experimental results indicate that the power 
measurement error due to sampling of sinusoidal voltages and 
currents becomes excessively large when the number of obser- 
vation periods approaches one-half the size of the microcomputer 
data memory allocated to the storage of either the input sinu- 
soidal voltage or current. 
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I. INTRODUCTION 


This report contains a study of power measurement errors 
due to the bandwidth and phase mismatch of a power meter and 
the sampling of the input voltage and current of a power 
measuring instrument [1,2]. 

The power measurement error due to the frequency response 
of a wattmeter is general to all types of wattmeters both 
analog and digital. The frequency response of the power 
measuring instrument is affected not only by its internal 
electrical system but also by the external electrical system 
that exists between the voltage and current that is to be moni- 
tored. For example, the frequency response and phase shift of 
the current shunt and voltage divider which is an integral part 
of the power measuring system should be considered when studying 
the power measurement error. Even electrical cables connecting 
the current shunt and voltage divider to the wattmeter should 
be considered as part of the power measuring system if the 
electrical cables are appreciable part of a wavelength. Not 
only do the cables introduce amplitude and phase distortion but 
can introduce time delays between the voltage and current. The 
exact behavior of the frequency response will be dependent on 
the choice of wattmeter, voltage divider, current shunt, end 
the electrical connections. 

In order to perform a power error analysis, the response 
of these electrical elements must be specified. In this 
investigation, it will be assumed that the transfer functions 
are first order between the points where the voltage and current 
are measured and where the multiplication of the voltage and 
current is actually performed. Hence the frequency response 
of the transfer functions decreases approximately at the rate 
of 20 db/decade beyond the corner frequency of the transfer 
function. 
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If the input voltage and current to the power meter are 
sampled, another power measurement error is introduced due to 
sampling. It will be shown that if the number of observed periods 
of a sinusoidal voltage or current is exactly equal to one-half 
the size of data memory allocated to the storage of either the 
input sinusoidal voltage or current, the power measurement error 
will be so large that the power measurement will be meaningless . 
This result is based on the assumption that both data memory 
size and the sampling rate are finite. Under the condition that 
the number of observed periods is one-half the data memory size, 
the number of samples per period is exactly two, which is exactly 
the minimum number of samples required according to Shannon's 
theorem [3] . 

The power measurement error due to the bandwidth and sampling 
cannot be analyzed until the current and voltage waveforms are 
specified as a function of time. In other words, the behavior 
of the power measurement error will change with the waveform of 
the input voltage and current. It can be easily shown if the 
input voltage and current are assumed periodic functions of time 
and superposition of harmonic power is applied, the total average 
power, P.J., can be represented by summing the harmonic average 
power, P n 


P 


T 


n ! 


* P 

0 n 


( 1 ) 


The normalized error of the total power is related to the nor- 
malized error of the n^ harmonic average power and assumes the 
form 


4 V P I- Jo <V P T>“V P n>- < 2 > 

Note that the normalized error of the n C ^ harmonic average power 
is multiplied by a weighting factor, P N /P,p. In order to determine 
this weighting factor, the waveform of the voltage and current 
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must be specified. In this study of power measurement error, 
two waveforms will be used; they are: the sinusoidal signal 
and the waveforms generated by a simplified DC/ DC chopper 
system. 

The power error analysis using a single frequency approach 
offers the advantage of studying the behavior of the power 
measurement error using trigonometric equations which can be 
manipulated into mathematical forms that are easily understood. 
The results of the single frequency error analysis can be 
extended to more than one frequency, but the error analysis 
becomes more complex. 

The waveforms generated by a simple DC/ DC chopper system 
offer the other extreme in error analysis when a complete 
Fourier analysis has to be performed in order to study the 
behavior of the power measurement error. The choice of simple 
DC/DC chopper system was based on the fact that such a system 
is used to control large amounts of power such as in an electric 
vehicle. Since electric vehicle's efficiency is very important 
from a practical viewpoint, an understanding of the power 
measurement error as a function of system parameters has a 
very high priority if the vehicle's driving range is to be 
optimized. 

This study also investigates the feasibility of using a 
microcomputer as a wattmeter. The results of the investigation 
show that if the number of periods observed by the microcomputer 
approaches half the data memory alloted to either voltage or 
current, the microcomputer's output has a large power measure- 
ment error. This is due to the fact that the microcomputer 
samples the input signals at a finite rate. 
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II. ELECTRICAL MODEL OF A PHYSICAL WATTMETER 

The power delivered to an electrical system as a function 
of time is given by the product of the instantaneous voltage 
and current at the input terminals to the electrical system. 
Since the voltage and current will in general vary from in- 
stant to instant, the product of instantaneous voltage and cur- 
rent will also be a function of time. This product is de- 
fined as instantaneous power, 

p(t) - e (t) i(t) (3) 

where e and i are the instantaneous voltage and current. 

The average value of the instantaneous power or total average 
power is defined as 


P 


T 


_1 

t 2 -ti 



P(t)dt 


(A) 


where t^ and t 2 represents the beginning and end of an 
observation time. Since no restriction is placed on either 
the instantaneous voltage or current, the average power will 
in general be dependent on the length of observation time. 

If the instantaneous voltage and current are periodic, 
the total average power is 


1 f P 

P T " TT J P (t)dt 


P 'o 


(5) 


where Tp is the period. For the general case where voltage 
and current are periodic, Equation (5) can be used to calcu- 
late total average power providing that Tp represents the 
time of the fundamental period. 

Employing Fourier series analysis to a periodic voltage 
and current, the complex waveforms can be expressed in terms 
of an average component and an infinite sum of cosine functions. 
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The general trigonmetric form for Che Fourier series [4] is 

OD 

e(t) - V + z V_cos (n<u rt t +ip_) (6-a) 

0 n-1 n on 

«Q 

i(t) - I c I n cos (nu 0 t +* n ) (6-b) 

where V Q and I Q represent the average value of e(t) and i(t) 
respectively; V , and I n> $ n represent the magnitude and 
phase angle of the n c ^ harmonic of voltage and current, respec- 
tively; and u) Q represents the fundamental angular frequency 
(« 0 - 2 tt/T p ). 

Using Parseval's theorem [5], it can be shown that the 
total average power for a general periodic signal can be ex- 
pressed as an infinite sum of harmonic average powers. The 
result is 

oo 

P T * Vo + % n^l Vn cos (7 > 

Equation (7) relates the harmonic average power content of a 
periodic signal to the total average power. Note that for 
the n c ^ harmonic both V R and I could be large in magnitude and 
still the average power associated with the n c ^ harmonic would 
be zero if ^ - $ n ■ +90°. This point will be emphasized in 

a later section. 

Although the above equations may describe an ideal power 
meter, physical wattmeters deviate from the ideal case because 
the power instrument will measure quantities that are pro- 
portional to input voltage and current. If the proportionality 
constants are known, they can be accounted for in the design 
of the wattmeter so that the instrument reads the correct value 
at any frequency. 

A more general approach to modeling a wattmeter is shown 
in Figure 1. The transfer functions from the respective voltage 
and current inputs to the multiplier account for any electrical 
networks that exist between the inputs of the power meter and 
multiplier. For example, a current shunt could be considered 
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a part of Hj and a voltage divider could be included in the 
description cf Hy. Besides the external networks such as the 
current shunt and voltage dividers, internal electrical net- 
works to the multiplier must also be considered as part of 
Hj and Hy. For purpose of analysis, it is assumed that the 
transfer functions describe first order systems that can be 
represented by 

Hy ■ 1/(1+ (Jn« 0 )/« 1 ) (8-a) 


H x - 1/(1 + (jknu. 0 )/ui 1 ) 


(8-b) 


where represents the angular frequency where the magnitude 
of Hy is 3db down from the mid-band gain. This angular fre- 
quency is sometime.® referred to as the 3db angular corner fre- 
quency or just corner frequency. The parameter k allows the 
two transfer functions to have different comer frequencies. 
The factor nu> 0 represents the set of discrete angular fre- 
quencies associated with nonsinusoidal periodic signals as 
determined from the Fourier series. 

Since the transfer functions have magnitudes and phase 
angles, both quantities will affect the measurement of the 
n fc k harmonic average power. Equation (7) can be modified to 
account for the affects of Hy and Hj on the measurement of 
total average power. Since H v and Hr affect the magnitude and 
the phase of the n n harmonic of voltage and current, re- 
spectively, the approximate total average power [6] which 
will be indicated by the power meter is 


APP 


Vo + * n-1 


£ V n I n cos U n - $ n +tan"^ (n<j 0 /aj^) -tan"* (kn^/w^) ) 


((l+(nu, 0 /a) 1 ) 2 ) (1+ (knu) 0 / w 1 ) 2 ) 


(9) 


Equation (9) shows that as long as is finite, there will 
be an error in the measurement of power because the results 
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predicted by Equation (9) will be different from the results 
predicted by Equation (7) . The magnitudes of the transfer 
functions affect the product of V n I n while the corresponding 
transfer function phase angles affect the resultant power 
phase angle. If k«l, the resultant power phase angle is in- 
dependent of u because the two arc tangent expressions in 
Equation (9) cancel each other. In this case, only the 
magnitude of the transfer functions introduce an error in the 
measurement of average power. As the bandwidth of the instru 
ment decreases (o>^ -*> u Q ) , the power measurement error will 
correspondingly increase until errors are so large that a 
power measurement will have little or no meaning. 

If the input signals to the wattmeter are periodic, 
the bandwidth of the wattmeter should be large enough to en- 
compass all the important power contributing harmonic com- 
ponents in order to minimize the power measurement error. 
Section III will develop the relationships between bandwidth 
of a wattmeter and the corresponding power measurement error 
for a specific set of periodic signals. 
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III. ERROR ANALYSIS IN THE MEASUREMENT OF AVERAGE 
POWER OF SINUSOIDAL SIGNALS 

This part of the report will primarily be devoted to 
the error analysis in the measurement of average power of 
sinusoidal signals. Although single frequency power error 
analysis is somewhat restrictive, it does offer a clear view 
about the power error behavior. The results can be extended 
to include the non- sinusoidal periodic case if the true 
harmonic power content is known. 

A. Power Measurement Error Due to Transfer Functions H„ and Hj 

As discussed in Section I, the two transfer functions Hy 
and Hj affect wattmeter accuracy especially when the wattmeter's 
bandwidth is of the same order as the highest important harmonic. 
The actual distribution of harmonic power will not only depend 
on the amplitudes of the harmonic voltage and current, but 
also on the phase relationship between the voltage and current. 
Thus, the distribution of harmonic power will depend strongly 
on the profile of instantaneous voltage and current in the 
time domain. Since there are an infinite number of voltage 
and current profiles, the error analysis does not lend itself 
to a closed form analysis unless specific profiles are chosen 
for the instantaneous voltage and current. A sinusoidal wave- 
form offers the simplest approach because the power equations 
can be expressed in closed form through the superposition of the 
power theorem, the results can be extended to more complex 
periodic signals, because Equation (9) is general for any 
periodic waveform, one term out of the infinite set of power 
terms can be used to represent the sinusoidal profile and is 
given by 

P APP " ^m 1 !^ 2 ^ CO- 4 * 2 ) d+(kx) 2 ))^) (cos(0+tan" 1 (x)-tan“ 1 (kx))) (10) 
where V m and I m represent the amplitude of the sinusoidal voltage 
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and current., respectively; x represents the normalized angular 
frequency (x - u/u^); and e represents the phase angle be twee;* 
the voltage and current. The parameter k has the same meaning 
as in Euqation (8). The true average power is given by the 
relationship 

P T - limit P APP » (V_I m /2)cose. (11) 

1 r m 141 

k finite 


If Equation (10) is divided by Equation (11) , the re- 
sultant function will be independent of the peak voltage and 
current and a funct ion of x,k, and e only. The ratio p aPP^ P T 
is defined as the normalized correction factor, NCF. The true 
average power P T is equal to P App/ NCF * Therefore, it is 
important that NCF be as near to unity as possible for 
accurate power measurements. The normalized correction factor 
NCF is given by 

NCF-(1/ ((1+x 2 ) (l+(kx) 2 ))^) (cos (e+tan'^x) -tan" 1 (kx)) /cose) . (12) 

In general, NCF will approach + • as the phase angle approaches 
+90° respectively, providing x is not zero. Hence the 
power measurement becomes meaningless for very small power 
factor? if tan _1 x - tan"“(kx) V 0. Equation (12) shows that 
if k is equal to unity, NCF will be independent of the phase 
angle. By carefully matching the frequency response of Hy 
and Hj, (k»l), Equation (12) will be a function of the normalized 
frequency x, only. 

Figures 2, 3, and 4 show the behavior of the correction 
factor as a function of the normalized frequency with k and 
e as parameters. Except for the case where k"l, the correction 
factor, NCF, is strongly dependent upon the phase angle and 
becomes more pronounced as the power factor approaches zero. 
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Accord ing to Equation (12), the correction factor should dis- 
play some sort of oscillatory behavior as a function of the 
normalized frequency, with oscillations becoming more pro- 
nounced at low power factors. The normalized frequency in 
Figures 2, 3, and 4 is restricted between zero and unity. 

Above unity the normalized frequency implies that the actual 
frequency exceeds the bandwidth of either Hy of H^; this is 
not the normal frequency range for the power meter. Although 
Figures 2 and 3 are informative, they are not very useful 
from a measurement viewpoint because the correction factor 
varies over too large a range for meaningful results . By re- 
stricting the range of the normalized frequency (see Figure 4) , 
a tighter correction factor tolerance (0.98 £ NCF <_ 1.02) can 
be achieved. This, of course, restricts the useful frequency 
range of Hy and Hj. 

It is obvious from Figure 2 that Hy and Hj should have 
identical frequency responses (k=l) if the correction factor 
is to be independent of the phase angle. For x -*• 0, NCF •+■ 1 
regardless of the values of k and 0. For k = 1 + Ak, which 
allows for slight differences in the frequency response be- 
tween Hy and Hj, a linear equation representing NCF as a function 
of x, assuming x-*-0 with e and Ak treated as parameters, can 
be derived indirectly from Equation (12) by expressing 
Equation (12) as a fourth order polynominal . The polynominal 
is given by 

(NCF)k 2 x 4 +( (NCF) (1+k 2 ) -k)x 2 +( (1-k) tane )x+NCF-l=0 . (13) 

For the case k = 1 + Ak, x + 0, NCF -*■ 1, Equation (13) is 
reduced to the linear equation 

NCF ■ l-(Aktane)x. (14) 

The slope, Ak tane, becomes very large in magnitude as the 
power factor angle approaches +90° for Ak40, and is positive 
for 0 negative and negative for & positive. The result is 
consistent with Figures 2, 3, and 4. 
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For the case where the voltage and current are periodic 
functions of time, the true total power can be expressed as 
a function of approximate harmonic power and the corresponding 
NCF n by 

P T ■ io P APP n /[)Cf n- < 15 > 

Although Equation (15) gives the same results as Equation (7) , 
the form of Equation (15) is such that each approximate n C ^ 
harmonic power must be known or measured and the corresponding 
n n narmonic phase angle be known or measured in order to 
determine NCF n for the corresponding n C ^ harmonic. However, 
Equation (15) does emphasize that a wattmeter having a small 
power measurement error must have NCF n as close to unity as 
possible with a correspondingly small normalized frequency. 
This can be accomplished using a wattmeter which has a band- 
width much broader than the important harmonics contained in 
the voltage and current signal and one which has identical 
frequency responses for Hy and (k=l) . If the bandwidth 
of the wattmeter is many times larger than the important 
harmonics being monitored by the wattmeter, the frequency 
matching of the two transfer functions, Hy and Hj can be re- 
laxed to some degree. 

B. Power Measurement Error Due to Sampling the Voltage and 
Current 

This portion of the report will primarily be concerned 
with power measurement error analysis associated with sampling 
voltage and current and the determination of total average 
power based on these samples. For example, this type of 
power measurement error occurs when a microcomputer is used 
as a power meter; the input analog signals are sampled, con- 
verted to digits, stored, and manipulated to generate an 
output equal to the average power. In order to develop a 
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power error equation in closed form, it will be assumed 
that the input voltage and current are sinusoidal. 

Before investigating the behavior of power measure- 
mei:t error due to sampling, it is instructive to first ex- 
plore the power measurement error associated with continuous 
signals. This would be equivalent to sampling at an infinite 
rate with a power computing system having an infinite memory. 
Instantaneous power can be expressed as 

p(t) » (V m I m /2) (cos 9 -cos(2ut +e )). (16) 

Using Equation (4) for the definition of average power and 
letting observation time interval, the approxi- 

mate average power due to the finite T Q is 

P * (V m I m /2)cose-((V m I m T p )/(B7rT o ))(sin(4TrT o /Tp+0)-sin0) (17) 

where Tp = 2ir/oi. The true average power is given by the first 
term in Equation (17) . The second term can be interpreted 
as an error term which oscillates with a period Tp/2 and de- 
creases as the observation time, T Q , becomes large with re- 
spect to the period Tp. The normalized power error, NPE, is 
defined in terms of true average power, 

NPE = (P-P)/P T = (sin(4irC+0) - sine)/ (4ttCcos9 ) (18) 

where C = ^ 0 ^P' t * ie num ber of periods or partial periods 
that are observed in time T Q . The normalized power error is 
a decreasing oscillatory function of period C and is zero 
for multiples of C/2. For a given C (C not equal to a 
multiple of %) , NPE approaches + « as 0 approaches +90°. 

Using L' Hospital's Rule, it can be shown that if C is equal 
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to a multiple of %, NPE approaches zero even as 9 approaches 
+90°. 

It can be concluded from the above discussion that a 
wattmeter capable of observing C periods equal to a multiple 
of % periods, the NPE will be zero for any value of 0. If 
C is other than a multiple of % periods, NPE is not zero and 
there will be a difference between the true average power 
and approximate average power. This power error is usually 
small when using a conventional laboratory wattmeter because 
the observation time is usually much greater than the period 
of the signal (human response time must be considered) being 
monitored by the wattmeter. However, a microcomputer, with 
a very large memory but finite sampling interval could be 
programmed to determine the approximate average power over an 
observation time interval comparable to the period of the 
input voltage and current. The value of NPE could become 
significant, especially if the phase angle approaches +90°. 

If rhe input voltage and current are sampled at a finite 
rate and the memory size is finite, the above results for the 
continuous case are no longer valid in a strict sense. The 
normalized power error no longer keeps on decreasing with 
increasing C because of the finite sampling rate and finite 
data memory size. This condition is observed when a micro- 
computer with finite data storage capability is programmed to 
function as a wattmeter. The analog input signals are sampled 
and processed to obtain the approximate average power. The 
gating time or observation time will depend on the memory size 
and time between samples and is given by 

T * NT (19) 

o s 

where N represents data memory size and T g the time between 
samples. For a given data memory size, the observation time 
is proportional to the sampling time interval. For large 
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sampling frequencies (T is small) it is possible for T 
to become comparable to the period of the input signals 
thus causing the power measurement error to be excessive. 

The discrete form of Equation (16) can be written 

p(nT g )»(V m I m /2) (cos 9 - cos (4irnT s /T p +0)). (20) 

The corresponding approximate average power is found by 
summing Equation (20) over memory size N. P is given by 

N 

¥ - (1/(N+1)) n l Q p(nT s ). (21) 

It is convenient to define a power measurement error simi- 
lar to Equation (18) except that it will be called the nor- 
malized sampled power error, NSPE. The result is 

N 

NSPE = 1-P/P = (1/ (N+l) ) ( n ^ Q cos(4TrnT s /T p +0)/cos0) . (22) 

Expressing the cosine function in terms of exponentials, 
and recognizing that the experiential summation is a geometric 
series, the normalized sampled power error can be expressed 
in term of periods or fractional periods C. The result is 

NSPE*(1/ (2 (N+l)cos 0 sin(2irC/N) ) ) (cos (0+*) -cos (0-H|»+4 itC) ) (23-a) 

* » tan* 1 (sin(4irC/N)/ (2sin 2 (2irC/N))) . (23-b) 

If N is allowed to approach infinity, 

If/ -*■ IT / 2 

cos(0+^) -*■ -sino 


cos(0+ +4irC)->- -sin(0+4irC) 
sin(2irC/N) 2irC/N . 
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The result is exactly the same as for the continuous case 
(Equation (18)). If C is chosen to be a multiple of %, NSPE 
is zero. This is also consistent with the continuous case. 

The major difference between Equations (23 -a) and (18) 
is the factor sin(2irC/N) instead of 4irC, The sine function 
is cyclic and appears in Equation (23-a) because the input 
signals are being sampled at a finite rate and the memory 
size is finite. For a specified N (memory size) and N >> 1, 
the NSPE will approach +1 and -1 as C approaches 0 and N/2, 
respectively. The case for C approaching 0 implies that the 
observation time as well as the sampling time is approaching 
zero (Equation (19)). This means P/P is approaching zero in 
the limit or NSPE is a +1 (Equation 22)). As C approaches 
N/2, the number of samples per period, N/C, approaches 2. 

The angle <f> (Equation 23-b) approaches tt/ 2 and NSPE (Equation 23-a] 
becomes indeterminant. Using L' Hospital's Rule, it can be 
shown that NSPE -*>-1 as C + N/2. 

If C in Equation (23) is replaced by C n + AC, where C R 
is an integer, the behavior of NSPE can be studied within 
a given period C n . It is assumed that the scanning parameter, 

AC, is restricted to a range from zero to unity. This pro- 
cedure permits the comparison of NSPE of one period with 
another period as shown in Figure 5 and 6. Figure 5 shows 
the behavior of NSPE for C n = 1 and C n * 255 assuming N * 512 
and 9 = 30°. The plot shows for C ■ 1 that the NSPE will 
reach a maximum of approximately 0.04 and a minimun of approxi- 
mately -0.10 for AC equal to approximately 0.10 and 0.35, 
respectively. If the power instrument observation time happens 
to coincide with C = 1.15, or C * 1.35, the power measurement 
error would be too large for meaningful measurements . For 
the period 255, which is just before the period where C equals 
512/2, the NSPE begins to increase because of the sin(27rC/N) 
factor in Equation (23-a). Note that NSPE for C = 255 is read 
along the right hand side of the vertical axis in Figure 5. 
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Again, the power measurement could be meaningless if the 
observation time is chosen incorrectly. In both cases, 
the phase angle was fixed at 30°. 

Between C n =* 0 and C n = 256, NSPE decreases because 
the sin(2irC/N) factor in Equation (23-a) is increasing in 
value. Figure 6 shows the behavior of NSPE for C n * 50, 

100, and 127; e - 30°; and N * 512. Note that the NSPE is 
considerably less than for C n - 1 and 255. According to 
Figure 6, if 50.4 periods were observed, the value of NSPE 
is approximately -0.0017 and the approximate average power 
is P t ( 1+.0017) or 0.17% above the true power. This power 
measurement error is well within the acceptable limit for 
most power measurements . Figure 6 also reflects that the 
maximum absolute error is smaller for C n = 100 rather than 
for C n » 127 where the sin(2*C/N) is approximately unity. 

This is due to the phase shift angle e which is 30° for this 
case. Since NSPE is oscillatory, it would be advantageous 
to determine the maximum and minimum values for any period 
between C=1 and 255. This approach would give the worst case 
error for any period (1 < C < 255) . The behavior of the 
minimum worst case as a function of C can be studied for 
various values of phase angles. 

A computer program was developed to determine NSPE ma y 
and NSPE m ^ n for any period C n where C n could vary from 1 to 255. 
Figure 7 shows the behavior of NSPE may as a function of periods 
observed for e = +30° and N * 512. For the stated phase angles, 
NSPE ma x reaches approximately the same global minimum but 
in different numbers of periods. This minimization is pri- 
marily due to the sin(27rC/N) factor and the phase angle. 

Figure 8 shows the variation of NSPE m ^ n for the same range of 
C as in Figure 7. 

The results of the section reflect the importance of having 
a matched frequency response for Hy and when the bandwidth 
of the power meter is comparable to the highest important 
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harmonic power associated with periodic voltage and current 
signals. Also, harmonics with large voltage and current 
magnitudes may contribute very little average power if the 
phase angle is near +90°. 

In summary, it was found that sampling the voltage and 
current and processing the results using a microcomputer as 
a wattmeter can generate excessive power measurement errors 
if the observation time is comparable to the period of the 
input signals or if the number of observation periods approaches 
half the memory size used to store the samples of voltage or 
current in the computer . 
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IV. POWER ERROR ANALYSIS OF A SIMPLE CHOPPER SYSTEM 

As pointed out in Section II, power measurement error 
is strongly affected by the waveform of the instantaneous 
voltage and current. Hence, in order to understand the be- 
havior of power measurement error, a specific set of waveforms 
must be investigated. This section of the report will analyze 
the power measurement error associated with an ideal simple 
chopper system [7] that controls the power from the source to 
a load. In this case, the load will be an ideal shunt motor. 

The duty cycle of the chopper and the motor's back emf have 
strong influences on the power measurement error. 

In order to evaluate the harmonic input and output power 
content in an electronic chopper or power controller such as 
employed in an electric vehicle, the steady state instantaneous 
input/output voltage and current must be known in terms of 
the system parameters. A simplified schematic for a DC/DC 
chopper is shown in Figure 9. Power that is delivered to 
the shunt motor from the DC source is controlled by switch 
SW that is activated by an electronic system not shown in the 
model. By varying the on-time with respect to the time of 
one period, power supplied to the motor can be varied from 
zero to a maximum value which depends on the system parameters. 

For the purpose of analysis, it will be assumed chat the time 
constants of the electric motor and the connecting mechanical 
load are much larger than the time constant of the power controller. 

The resistances R^, R 2 . and R^ represent the internal 
resistance of the DC source and input leads; on-resistance of 
the switch SW; and resistance of the armature of the shunt 
motor, respectively. Inductance accounts for the inductance 
of the armature. Voltages E^ and E 2 represents the open-circuit 
voltage of the D. C. source and the back emf of the motor, 
respectively, while voltages E^ and E^ account for the semi- 
conductor offset voltage. Listed in Figure 9 are the numerical 
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values of the system parameters used In the harmonic analysis 
of the power controller's input and output signals. 

Diode is necessary to allow the current i£ to con- 
tinue to flow after SW is turned off; otherwise the induced 
voltage across will approach infinity as i 2 decreases 
instantaneously to zero. This will cause a voltage breakdown 
in the power controller. The free-wheeling action of the 
diode and increases DC component substantially by taking 
advantage of the energy stored in the inductor. 

Typical voltage and current waveforms are shown in 
Figure 10. The on-time is represented by Tq and the period 
by Tp. Obviously, T^ can be varied from 0 to Tp and for 
purpose of analysis, the period is fixed at 0.001 second. 

The harmonic structure of the input/output voltage and 
current will be required in order to study the behavior 
of the harmonic power content. The on-time and back emf are 
the two system parameters that can be varied easily In an 
electric vehicle. The duty cycle, T^/Tp, is controlled by 
the vehicle's operator via the power controller while the 
back emf depends on the magnetic flux and the speed of the 
motor. Neglecting armature reaction, the net flux will 
essentially remain constant for a given shunt field current 
and armature current. Under this condition, the back emf is 
directly proportional to the speed of the armature. 

The amplitude and phase angle of the harmonics contained 
in e^, i^, and ^ will vary with duty cycle for a given 
back emf. Before quantitatively expressing the harmonic 
behavior as a function of the on-time, a Qualitative assessment 
can be conducted. This approach reinforces the quantitative 
results. Starting with a 1007. duty cycle, the voltages and 
currents in the system are pure DC assuming the electrical 
system in Figure 9 nas reached a steady-state condition. 

As the duty cycle is reduced from 1007., the signals become 
rich in harmonics and this richness grows as the duty cycle 
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decreases. However, as Che duty cycle approaches zero, Che instan- 
caneous volCages and currencs and corresponding harmonics as 
well as Che DC componenc approach zero. For a given ducy 
cycle, Che harmonic concenc in Che inscancaneous volcages and 
currencs decreases as Che back emf approaches che input vol- 
Cage. This is obvious because Che back emf bucks che volcage 
E^. The inscancaneous volcages and currencs are cherefore very 
rich in harmonics when che ducy cycle is small and che back emf 
is approaching zero. 

Since Che power associaCed with any given harmonic de- 
pends on che amplicude of che volcage and currenC as well as 
che cosine difference of che phase angles, large harmonic 
ampliCudes ChaC are approximacely in quadracure wich each ocher 
concribuce very liccle Co Che CoCal average power. If che DC 
componenCs of che inscancaneous volcage ::nd currencs are sup- 
pressed as shown in Figure 11, very inCeresCing resulcs can 
be observed. The fundamencal inpuc volcage and currenC as 
shown by Che docCed curves are approximately 180° ouc of 
phase. The exacc phase relacionship will be a funcCion of che 
cuCy cycle, However, Che fundamencal ouCpuC volcage and cur- 
renC appear Co be in quadracure, Chus contribut, '<? very liccle 
Co Che CoCal average power. This is primarily ■ Co Che in- 
ducCive filCering acCion by Che armacure inductance which cakes 
place via Che diode D^. QualiCacively , che frequency specif i- 
caCions of a waCCmeCer monicoring Che inpuc power of che 
power conCroller will be more scringenc chan for a wattmeter 
Chat is monicoring che power ouCpuC of che same conCroller. 

A. Fourier Analysis of Che Signals Associated with che Simple 
Chopper SysCem 

In order Co evaluace Che harmonic concenc of e^, i^, e£ 
and ± 2 > Che syscern muse be in steady state; otherwise the 
Fourier series analysis Is meaningless. A set of equations 
representing the instantaneous, voltages and currents will be 
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derived based on circuit boundary conditions. The results of 
such analysis will allow the determination of total average 
power, harmonic power content via the Fourier series, and 
power error due to a wattmeter with a finite bandwidth. 

The current, i^, can be expressed in the form 

11 ■ Ki + l^expf.-t/T^) ; 0<.t <_T D (24-a) 

i 0 - ° ; T d < t<T p (24-b) 

where represents \he value of i^ if the switch SW remains 
closed indefinitely and is equal to (E^-E^-E 2 >/ (R^+R^+I^ * 
K 2 " i 2 ^”^” K l’ an< * T^/ (R 1 + R 2 +R 3 ) . The current ^(0") is 
the minimum value of the output current as indicated in 
Figure 10. During the time interval when SW is closed, i^*^, 
and during the remaining portion of the period when SW is 
open, i 2 can be expressed as 

1 2 - K 3 exp(-(t-T D )/T 2 ) ; T D < t<T p (25) 

where ^ an< * *- s assume d that i 2 never 

reaches zero beFore the switch is closed again. This assumption 
is reasonable since the time constant T 2 is usually much 
larger than the period T p . Equations (24-a) and (25) must 
have the same value at t«T D and Equation (25) must be equal 
to ^(O") at T«T p since the system is assumed to be in a 
steady state condition. Thus, 

K x + (i 2 (0*)-K 1 )exp(-T D /T 1 ) - K 3 ; t-T D (26-a) 

i 2 (0") - K 3 exp(-(T p -T D )/T 2 ) ; t-T p (26-b) 

From Equations (26-a) and (26-b), ^(O”) can be expressed 
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in terms of Kj_, Tp, T p , and T 2 as 


i 2 (0‘)-K 1 (l-axp(-T D /T 1 ))/(exp«T p -T D )/T^exp(-T D /T 1 )). (27) 


The vs.lue of i 2 (0") at C *G and t«T p is in agreement with the 
electrical bounds on the system, namely, i 2 ( 0")"0 and 
^(O - ) ■ ■ (E^-E^-E 2 >/(R^+R 2 + R 3 ) . respectively. The instan- 

taneous input and output voltages can be expressed in terms of 
the respective currents i^ and ^ 


e l - E 1 - R 1 i 1 


e l “ E 1 


; 0 < t < T d 


: T D r 1 T P 


®2 m “ (R^+R2>i^ ; 0 t <_ Tj 


e 2 " ” E 3 


•- T D < t < T p 


(28-a) 

(28-b) 

(28-c) 

(28-d) 


Equations (25), (26), (27), (28) are the necessary eet of 
equations to determine the total input and output power, 
voltage and current harmonics, harmonic power content, and 
power error due to a wattmeter with a finite bandwidth. 

The harmonic content associated with the controller's 
instantaneous input and output signals can be expressed as 
a set of algebraic equations whose coefficients are a function 
of the system parameters. The equations are listed below. 

See the Appendix for the functional relationship between X^, X 2 , 
Yp Y 2 , Y^,Yj and the system parameters. 

Input Side of the Power Controller 


Voltage Terms 

W o" E l T D /T P -R l K l T D /T p- (R l T l K 2/ T P )(1 ‘ ex P ( ’V T l ))+E l (T P’ T D )/T P (29-a) 
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W x - -RCK^ + K 2 X 2 ) (29-b) 

W 2 - -^(K^ + K 2 Y 2 ) (29-c) 

Current Terms 

W oo " K 1 T D /T P + (^/Tp) (1-exp (-Tp/Tj^)) (30-a) 

W 3 - + K 2 X 2 (30-b) 

W 4 - + K 2 Y 2 (30-c) 

|V I " << w i> 2 -K W 2> 2 > % and - -can" 1 (W 2 /W 1 ) <31-a) 


1 1 n 1 l - ((W 3 ) 2 +(W 4 ) 2 )^ »nd/l ni .. tin -l (vw ^ (31 . b) 

Output Side of the Power Controller 
Voltage Terms 

X 0 K E i- E 4 )V T P“^i'«2^ K l T D /T P* ((R l' W 2 )T l K 2 /T P )(1 - e,c P ( - 1 b T i >) " E 3 (T P" T D ) ^ r P (32 ’ a 

Z 1 " ( e i" e 4) x i " (R 1 +R 2 )(K 1 X 1 +K 2 X 2 ) + E 3 X L (32-b) 

Z 2 - (E 1 -E 4 )Y 1 - (R 1 +R 2 )(K 1 Y 1 +K 2 Y 2 ) + E 3 Y 1 (32-c) 

Current Terms 

X bo" K l T D /T P +<K 2 T l /T P ) (l-expCT^Tj) )+(K 3 T 2 /T p ) (1-exp ( - (Tp-ty /T 2 ) ) (33-a 
Z 3 * K 1 X 1 + K 2 X 2 + K 3 Y 4 


(33-b) 
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Z 


4 



K 1 Y 1 + K 2 Y 2 + K 3 Y 5 

(33-c) 

- ((Z 1 ) 2 +(Z 2 ) 2 ) % and // / V no -tan" 1 (Z 2 /Z 1 ) 

(34-a) 

-(<Z 3 ) 2 +(Z 4 ) 2 )* and 

(34-b) 


Equation ( 31 ) and (34) were evaluated using a computer. A 
typical computer printout is shown in Figure 12. The first 
thirty harmonics of input/output voltage and current were 
calculated as a function of on-time T D for a given E 2 . For 
example, the case where T D * 0.0001 second and E 2 =25 volts, the 
input fundamental voltage and current are for all practical 
purposes 180° out of phase while the output fundamental vol- 
tage and current are 90° out of phase. This is consistent 
with the previous discussion about the harmonic phase relation- 
ship between voltage and current at both the input and output 
of the power controller. Although the value of the phase 
angles for the individual voltage and current are dependent 
on T D and E 2 , the angular difference between either the in- 
put voltage and current or output voltage and current remains 
independent of T^ and E 2 , which is 180° and 90° for the input 
and output, respectively. Results indicate that the amplitude 
of the harmonics decreases for increasing T^ and E 2 which is 
consistent with a physical system. The results of this 
section are needed for the calculation of the harmonic power 
content at the input and output of the controller. 

B. Total and Harmonic Power of a Simple Chopper System 

Total average input and output power can be determined 
from Equations (24), (25), (26), (27), and (28) by multipling 
e l ^ ®2 ky * 2 ' anc * averaging the results over one 

period. The results are 
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► 


P i -E 1 I 0i -(R 1 K 1 Z T D )/T p -(2K 1 K 2 T 1 R 1 /T p )(l-exp(-T D /T 1 )) 
- (K 2 2 T 1 R 1 /(2T p )Xl-exp(-2T D /T 1 )) 


(35) 


P o" <E r E 4 )I oi“< R l +r ^< K l 2T D /T P +(2K l K 2 T l /T P )(1 " ex P ( ’ T D /T l )) 
+ (K 2 2 T 1 /(2T p ))(l-exp(-2D/T 1 )))-(E3K3T 2 /Tp) <l-exp((Tp-T D )/T 2 )) 


(36) 


where I Q ^ represents the DC input current. See the Appendix. 

Harmonic distribution of average power can be determined 
from the knowledge of the magnitude and phase angle of voltage 
and current. Frequency roll-off due to the first order input 
•zransfer functions of a wattmeter can be incorporated very 
easily into the power equations. The results are 


P i 


APP 


W o W oo +(%) 1 £ 1 ^n^n^cos (-tan“^ (V^/W^)+tan”^ (W^Alj-t3n“''(nw ( 


N- 


“1 


+tan‘ 1 (kna, 0 /a) 1 ))/((l+(na> 0 /a) 1 ) 2 ) (l+(kna) 0 /a) 1 ) 2 ))^) (37) 

N 1 

P o “ Voo*<« E , CV no I no cos(-tan' 1 <Z 2 /Z 1 >4-tan" 1 (Z 4 /Z3)-tan" 1 (na» 0 /a, 1 ) 
APP n-1 

+tan" 1 (kn (1 j 0 /a) 1 ))/((l+(naj 0 /a> 1 ) 2 ) (l+fcna^/u^) 2 ))*) (38) 

where represents the highest harmonic of interest. The 
other terms were defined in Equations 6, 7, 8, and 9. As 
long as and k remain finite, the power calculated from 
Equations (37) and (38) will be an approximation of the 
results from Equations (35) and (36) even if •+■ «. This 
difference is due to suppression of the higher harmonic power 
terms by the transfer functions Hy and Hj. 

A computer program was developed to evaluate Equations 
(35), (36), (37) and (38) as functions of the back emf E 2 
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and the on-time T^. The computer program calculated a 
running power sum of approximate input /output power 
(Equations (37) and (38) as a function of N^, total 
input /output power Equations (35) and (36), and the power 
measurement error. 

If the transfer functions' comer frequencies are at 
least 100 times larger than the fundamental frequency of 
the periodic signals, it can be reasonably assumed that 
Hy and Hj are essentially unity and do not effect the power 
measurement. Assuming wide bandwidth transfer functions, 
the power measurement error will approach zero as the number 
of harmonics approaches infinity. However, from a practical 
viewpoint it is better to consider the necessary number of 
harmonics in order to achieve a certain power measurement 
error. For analysis purposes, a power error of 0.01 (17.) 
was chosen. With this power error objective, the desired 
error can be attained with relatively few harmonics assuming 
infinite frequency response. As the corner frequency of 
the transfer functions begins to approach the fundamental 
frequency of the periodic signals, more harmonics will be 
required to achieve a power error of 0.01. If the comer 
frequencies of the transfer functions are set equal to the 
fundamental frequency, the power error will be so large that 
it may be impossible to achieve a power error of 0.01 for 
certain combinations of back emf E£ and on-time t d- 

The results of power measurement error versus number of 
harmonics are shown in Figure 13 a-d. The figures are 
arranged in an array to emphasize the behavior of power 
measurement error and the required number of harmonics 
necessary to achieve a specified numerical value of the 
power measurement error as a function of E 2 and T^. The 
parameters k and f^ which set the comer frequencies of 
the transfer functions are varied to study the behavior 
of the power error. The fundamental frequency is 1000 hertz 
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(f Q * 1/Tp) . If an error of 0.01 cannot be achieved, this 
is reflected by indicating the actual error when all 30 
harmonics are considered. This will occur when f^ 

(f 1 = u ^/27r) is comparable to the fundamental frequency of 
the signal in the power controller. When Tp is equal to Tp, 
the power measurement error is zero because the power con- 
troller's switch is closed continuously and the voltages 
and currents are pure DC signals. This is shown in the last 
row of the array. 

The input power measurement error for the case of a 

wideband wattmeter is shown in Figure 13-a. With k=l and 

f,=100f , it is obvious that am error 0.01 can be achieved 
1 o 

for all combinations of E£ and Tp shown. It should be noted 
that the number of harmonics required to attain 0.01 error 
varies depending on the combination of E£ and Tp. The most 
severe case occuring when E2=25 and Tp=0.0001 second. 

Moving diagonally across the array, the required number of 
harmonics decreases for 0.01 error. This is consistent with 
the electrical model in Figure 9 since the duty cycle of 
the power controller is increasing and the E 2 is approaching 
E^ which is equal to 100 volts. 

Reducing the bandwidth of the wattmeter will cause the 
power measurement error to increase as shown in Figure 13-b. 

The first column depicts a condition where an error objective 
could not be attained even when thirty harmonics were considered. 
The result is due to the frequency roll-off of the transfer 
functions interacting with the higher harmonics. Again, moving 
diagonally across the array the 0.01 error can be achieved 
with less than 30 harmonics. 

If the bandwidth of the transfer functions equal the 
fundamental frequency, the power measurement error increases 
drastically and exceeds the 0.01 criterion over more than 
507c of the array as shown in Figure 13-c. It is obvious 
that a wattmeter whose bandwidth is equal to the fundamental 
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frequency would be practically useless, especially for low 
values of E£. 

The results shown in Figure 13-d are somewhat similar 
to Figure 13 -c except that the comer frequency of the two 
wattmeter transfer functions are numerically di-ferent. With 
k=*0.5, one corner frequency is twice the other. Power 
measurement errors are generally lower for k=0.5 which is to 
be expected since the effect of the transfer functions on 
the magnitude of the harmonics is less pronounced. However, 
with k not equal to unity, the phase angles of the transfer 
functions do not nullify each other according to Equations (37) 
and (38) and this causes a phase angle error which in turn 
affects the power measurement error. 

The results of this section of the report illustrate 
that the harmonic power content at the input is more important 
than that at the power controller's output. The reason is 
due to the electrical time constant of the shunt motor. In 
other words, the inductive reactance is the predominant con- 

pared to the armature resistance, R^, a^ the harmonic fre- 
quencies. Hence, a wattmeter having a broader bandwidth must 
monitor the power input to the power controller. 

Depending on the back emf of the shunt motor and the 
bandwidth of the power meter, there were particular cases where 
a power measurement error of 0.01 was exceeded, even with all 
30 harmonics. This occurred for low duty cycles and small 
back emf s . 
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V. POWER ERROR ANALYSIS OF A MICROCOMPUTER WATTMETER 


In Section II-B a theoretical investigation was per- 
formed to determine the power measurement error due to 
sampling an input sinusoidal voltage and current. The 
results indicate that the Normalized Sampled Power Error, 

NSPE, was very large if a portion of a period was sampled 
or if the number of periods sampled approached half of the 
data memory capacity, N/2, devoted to storing the voltage 
or current data. A global minimum value for NSPE occurred 
between C = 0 and N/2 with the position of the minimum NSPE 
functionally dependent on the phase angle 0. 

This section is primarily concerned with programming a 
microcomputer as a wattmeter and studying the behavior of the 
microcomputer's output error as a function of observed periods 
and phase angles. A commercial 16-bit microcomputer with a 
12 -bit I/O subsystem was programmed to accept voltage data on 
two of the 16 differential protected input channels. Input 
channels 0 and 1 were chosen in this case. After computations 
were completed by the microcomputer, the output channel 1 
displayed the results on a digital voltmeter. 

In order to be consistent with Section II-B, the micro- 
computer's data memory capacity was set at 512 for each chan- 
nel with channel 0 acting as the triggering channel. The 
triggering was accomplished by having the microcomputer check 
channel 0 for positive zero crossing. This permitted the 
microcomputer to process the equivalent of the same data in 
case the output reading was questionable. 

\ program flow chart describing the microcomputer 
operation as a wattmeter is shown in Figures 14 a-e. There 
are essentially three main parts to this program. The first 
segment of the program requires the programmer to set the 
sampling interval using a five digit number between 00010 
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and 10000. The microcomputer's internal timing system is such 
that a five digit number (xxxxx) divided by 46,875 equals 
sampling interval in seconds. Therefore, the sampling in- 
terval can be varied from approximately 213 microseconds to 
213 milliseconds. The minimum sampling time interval was 
chosen to be ten times the time between sampling channels 0 
and 1. 

Once the sampling time interval has been set and the 
teletype return carriage key actuated, the microcomputer 
observes channel 0 and waits until the signal goes through 
a positive zero crossing. At that point, both input channels 
0 and 1 alternately sample their respective input voltage 
signals, continuing until the data memories are completely 
filled. In this case, it would be 512 data points for each 
channel . 

The second part of the program multiplies the first 
datum from each channel and continues this multiplication in 
a sequential manner. After each multiplication, the result 
is added to the results of the previous multiplications with 
this operation continuing until there are 512 terms . The sum 
is averaged over the number products and scaled to fit the 
range of the analog output. 

The final part of the program directs the microcomputer 
back to the start position, ready for the calculation for the 
next sample interval. 

A mathematical relationship between the periods observed, 
C, and the input frequency, f, can be derived if the sampling 
time interval T g is specified. The result is 

C - NfT s (39) 

where N represents the size of the data memory. As previously 
stated in this section, the sampling time interval T g is 

T - (xxxxx) /46875 (40) 

s 
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where 10 <_ xxxxx <_ 10000 

Combining Equations (39) and (40), the number of periods 
observed is 


C - Nf (xxxxx) /46875 (41) 

A program following the flow chart in Figures 14 a-e 
was written and implemented on a commercial microcomputer. 

With N equal to 512 and sampling time interval capable 
of being varied from 213 microseconds to 213 milliseconds, 
the input frequency was set 36.74 hertz. The choice of 
frequency was based on the fact that small number of periods 
were to be observed in order to study the behavior pattern 
of NSPE. This low frequency also permitted the observation 
of NSPE for periods approaching N/2. 

The results of this section of the report are divided 
into two segments. The first segment deals with the compari- 
son of the calculated NSPE and theoretical NSPE for the same 
period of observation. The second segment studies the 
variation of the microcomputer's output as C approaches N/2. 

Figures 15 a-e compare the calculated NSPE with the 
theoretical NSPE for different phase angles. With an input 
voltage of 5.010 volts for both channels 0 and 1 and a 
frequency of 36.74 hertz, the correlation between calculated 
and theoretical NSPE is poor. This lack of correlation between 
the two normalized sample power error is believed to be due 
to either a slight input frequency shifting during the obser- 
vation time or jitter in the sampling time interval. These 
effects can cause an uncertaintly in the number of periods 
that are actually observed. The calculation of the theoretical 
NSPE is based on the numbers in the first column. The 
calculated NSPE tends to increase with increasing power factor 
angle which is consistent with the theoretical calculations 
of NSPE. 
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Figure 16 shows the microcomputer's output for different 
observation periods and phase angles . The bottom row marked 
true power output reading is used for comparison with actual 
output readings. For periods 265.64 and 269.66 the output 
readings are reasonably close to the corresponding theoretical 
output readings. However, at and near period 267.65, the 
output readings deviate substantially from the theoretical 
predictions. The differences are so great that the micro- 
computer would be useless as a wattmeter for periods observed 
in the vicinity of 267.65. It is noted that there is a 
slight difference between the period where the errors are 
excessive and the predicted period of 256 which is half of 
the data memory size (N«512) . This discrepancy amounts to 
approximately 57. and is believed to be caused by not knowing 
the absolute value of the input frequency accurately. As 
a cross check, the behavior of the microcomputer's output 
was studied in the vicinity two times 267.65 periods and it 
was experimentally observed at 536.5 periods the output error 
was excessive; this is consistent with theory. The ratio 
of 536.5 to 267.65 is 2.007 which is as close to two as 
experimentally possible. 

The results of this section demonstrate that a micro- 
computer can be used as a wattmeter providing that the periods 
observed for the computation is larger than several periods 
but less than half of the data memory size. As the periods 
observed approaches N/2, the error is so large that the out- 
put reading of the microcomputer is meaningless. This is 
consistent with the theory presented in Section II-B of 
this report where Equation (23) relates the functional be- 
havior of NSPE and the observed number of periods. In the 
denominator of Equation (23-a) the sine function, sin(2iTC/N), 
approaches zero as C approaches N/2 causing the power measure- 
ment error to be excessive. 
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Although the above analysis is based on a single 
frequency, the results are applicable for the general case 
of periodic signals. If superposition of power is assumed, 
it is possible for the observation periods of the harmonics 
to be near N/2 or multiples of N/2 causing a serious power 
measurement error. 
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VI. CONCLUDING REMARKS 

This study establishes the behavior of power measure- 
ment errors due to the bandwidth of H^ju) and the 

input transfer functions of a power meter, and the sampling 
of the input voltage and current of a power instrument. 

The frequency response of the power meter's input 
transfer functions affects both the magnitude and phase angle 
of the corresponding input sinusoidal voltage and current to 
the instrument. Since average power depends on the product 
of the magnitude of voltage and current times the cosine of 
the power phase angle, the effect of the frequency response 
is to cause both amplitude and phase distortion which affects 
the magnitude of the input sinusoidal voltage and current, 
and the power factor angle, respectively. The power measure- 
ment error due to H^ju) and HjCju) becomes more pronounced 
as the frequency of the input voltage and current approaches 
the bandwidth of H v (ju) and Hj.(ju>), and the power factor 
approaches zero. Assuming that H v (jw) and HjCju) are first 
order transfer functions, the power measurement error due 
to phase distortion can be eliminated if the two transfer 
functions have identical frequency responses . Reducing the 
power measurement error by matching the frequency response of 
HyCjuj) and HjCjtu) can be extended to higher order transfer 
functions. The reason this is true is due to the fact that 
the voltage and current phase error are subtracted when 
determining the actual phase angle. Of course, there will 
always be a power measurement error due to the amplitude dis- 
tortion of H v (ju) and H^Cju), and the behavior of this error 
will be strongly dependent on the order of the transfer 
functions . 

When the input sinusoidal voltage and current are 
sampled, a power measurement error can occur if the samples 
of the input voltage and current are taken over fractions of 
a period (not equal to a multiple of half periods) or if 
the total number of periods observed is equal to one-half 
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the memory size of either the voltage or current data memory. 
The behavior of power measurement error with respect to 
cycles observed is due, in part, to the cyclic behavior of 
the function, sin(2irC/N), which occurs in the demoninator 
of the Normalized Sampled Power Error, NSPE (Equation 23-a). 

As C N/2, the sine function approaches zero causing the 
power measurement error to increase. The numerical value of 
NSPE is bounded and equal to -1 as C •♦N/2 because the 
difference between the cosine functions is also approaching 
zero causing Equation (23-a) to have an indeterminate form 
at N/2. It can be shown that the bound of NSPE at C * N/2 is 
indeed -1 if N >> 1. 

Although the above power measurement errors are based on 
single frequency analysis, the results are applicable to more 
complex periodic signals. Unfortunately, in order to express 
the results in a closed form, the periodic signals must 
first be analyzed using the Fourier series before the har- 
monic power content can be determined. Once the harmonic 
power analysis is completed, each harmonic power measurement 
error can be determined based on the single frequency analysis 
with the resultant power measurement error being the sum of 
the harmonic power measurement errors. For periodic signals 
it would be more expedient to determine the harmonic power 
content contained in periodic signals by directly applying 
Fourier series analysis and using Parseval's theorem to 
determine the harmonic power distribution. This approach 
permits the inclusion of amplitude and phase distortion due 
to the presence of the transfer functions H v (ju>) and H I (jw). 

Since the total power measurement error is dependent 
on the weighting factor, P /P T , as indicated in the intro- 
duction of this report, the harmonic power distribution does 
influence the total power measurement error when the waveform 
of either the input voltage or current is changed. Therefore 
a power measurement error analysis must be performed on a 


» 
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given set of volcage and current signals. Input and out- 
put voltage and current signals of a simple DC/DC chopper 
system offer one such set of signals for power measurement 
error analysis. 

Using the duty cycle of the DC/DC chopper and back emf 
of a shunt motor which is the electrical load for the chopper 
as the independent variables, a power measurement error 
analysis was performed by computer. The results indicate 
that the worst case occurs when the duty cycle of the DC/DC 
chopper and the back emf of the shunt motor are numerically 
small. As a general rule, the power measurement error 
decreases as the duty cycle approaches unity and the back emf 
of the shunt motor approaches the DC/DC chopper supply voltage. 

The results also shows that if the shunt motor armature 
-'nductance is large with respect to the armature resistance, 
the output harmonic power content, excluding the DC power 
component, is small as compared to the input harmonic power 
content. This is due to the approximate quadrature phase 
relationship between the n c ^ harmonic output voltage and 
current which is caused by the dominance of the inductive 
reactance at the n c ^ harmonic. Hence, for the DC/DC chopper 
that was analyzed in this report, the bandwidth specifications 
for the wattmeter that is measuring the power input to the 
chopper would be more stringent as compared to the wattmeter 
that is measuring the power output of the chopper. Although 
not considered in the electrical model in this report, such 
electrical phenomenona as hysteresis and eddy-current 
losses can be represented by an equivalent resistance in 
series with the armature resistance. This resistance combi- 
nation may be sufficiently large as compared to the armature 
inductance to cause the output harmonic power to become 
significant. 



-38- 


This study also Investigated the feasibility of 
using a microcomputer as a wattmeter. Since the micro- 
computer samples the input voltage and current, the 
average power for sinusoidal input signals as determined 
by the microcomputer showed a power measurement error due 
to sampling. The power measurement error became significant 
when the number of observed periods C was approximately 
equal to N/2. The normalized sampled power error was 
substantially smaller when C was much different from N/2. 

The advantages of using a microcomputer as a watt- 
meter is that instantaneous sampled power, peek power, and 
the instantaneous sampled voltage and current can be stored and 
retrieved at a convenient time. When not using the micro- 
computer as a wattmeter, the microcomputer can be assigned 
other measuring tasks. 

The results of this study can be extended in several 
areas. For example, in the model of the simple DC/ DC 
chopper system, switching was assumed to be ideal with no 
losses during the opening and closing of switch SW. Switching 
transients tend to cause the rate of change of signals to 
be finite, thus decreasing the amplitudes of the high fre- 
quency components in the signal. However, if power analysis 
of a chopper is desired, the switching losses should be con- 
sidered. This is especially true if the fundamental chopper 
frequency is increased to a relative high value If the 
signals during switching can be represented by straight-line 
analysis, the switching power loss can be analyzed separately 
from the power analysis of the model of the DC/ DC chopper in 
this investigation. By combining the results of both power 
analyses, the behavior of the power measurement error would 
represent a more realistic chopper system. 

The ideal shunt motor that was assumed in this study 
could be modified to account for armature reaction and 
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eddy-current and hysteresis losses. The effect would be 
a larger output harmonic power content from the chopper 
as compared to the results of this Investigation. 

The behavior of the power measurement error with re- 
spect: to back emf and duty cycle would be of great Interest 
if the shunt motor were replaced by a series motor. The 
results could be useful in the electrical efficiency study 
of an electrical vehicle since most electric vehicles use 
a series motor because of its torque characteristic. How- 
ever, powe; analysis would involve a set of nonlinear 
equations due to the nonlinear relationship between the 
back emf of the series motor and the DC current output 
from the chopper. 

Although the use of a microcomputer in this study was 
restricted to positive zero crossing input signals, a com- 
puter program can be devised so that the microcomputer can 
accept a more general class of input signals. With such 
an arrangement, the output power reading could be auto- 
matically up-dated. Also, the microcomputer could store 
such quantities as peak voltage, current, and power and 
display this information with the appropriate command. 
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APPENDIX 


Fourier Analysis of the Input/Output Voltage and Current of 
a Simple DC/DC Chopper System and the Corresponding Input/ 
Output Total Average Power. 


In Section III of this report, a set of algebraic 
equations Equations (31) and (34) were presented based on 
the Fourier analysis of the instantaneous voltage and current 
at the input and output of a simplified DC/DC chopper s/stem. 
These equations were derived from the trigonometric form of 
the Fourier series and the general form is given by 


f(t) 


+ Z 
n=l 


n 


cos(na) Q t + l' n ) 


V 2 

a = (1/Tp) / f(t)dt 

° P -Tp/2 

Tp/2 

t a n = (2/Tp) J f (t)cosn<D Q tdt 

-Tp/2 

T p /2 

b n = ^ 2 / T p^ f f (t)sin naj Q tdt 
-T p /2 


(A 1-a) 


(A 1-b) 


(A 1-c) 


(A 1-d) 


n 


r 


n 


2 -1 
+ b and V = -tan (b / a ) 
n n n n 


(A 1-e) 


where f(t) is a periodic function; a Q represents the average 
component of f(t); a R and b R represent the in-phase and 
quadrature component of the n tlT harmonic, respectively; 

A n and ¥ n represents the amplitdue and phase angle of the 
n c h harmonic, respectively. The angular frequency, ui Q , is 
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equal to 2ir/Tp where Tp is the time of one period. 

This appendix will develop the Fourier series for the 
input and output voltage and current of the DC/DC simple 
chopper system and the corresponding input/output total 
average power. For convenience, the chopper ' s output signals will be 
analyzed first using Equation (A-l) and the representation 
of the instantaneous output voltage and current as 
developed in Section III. 

Output Voltage e 2 

~ E l” R l i l > 0 1 t <_ T D (A 2-a) 

®2 ~ ~ E 3 > i t i (A 2-b) 

i 2 * K i +K 2 ex P( _t / T i) ; 0 <_ t <_ T D (A 2-c) 


i 2 - K 3 exp (-(t-T D )/T 2 ) ; T D < t T p 


(A 2-d) 


a n *(2/T p ) J (E 1 -E 4 -(R 1 +R 2 ) (K 1 +K 2 exp(-t/T 1 )))cos(na Jo t)dt 

0 T p 

-(2E 3 /T p ) cos(nu) o T)dt 


(A 3-a) 



r Tp 

-(2E 3 /T p )J sin(nu) 0 t)dt 
T D 


(A 3-b) 


D “ 

a 0 -(l/T p )J (E 1 -E 4 -(R 1 +R 2 )(K 1 +K 2 exp(-t/T 1 )))dt-(E 3 /T p )/ dt. (A 3-c) 

0 T D 
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It is convenient to evaluate the following set of integrals 
because they occur in several places in Equation (A3) . 

T 

X 1 * < 2 / T p) jf P cos(nw 0 t)dt = 2sin(2imT D /T p )/(2Ttn) (A 4-a) 

T 

X 2 * (2/Tp) f D exp(-t/T 1 )cos(n< 1 . 0 t)dt - (2/ ((l/T^+O^) 2 )) 

° (A 4-b) 

((exp(-T D /T 1 ))((-l/(T p T 1 ))cos(2irnT D /Tp)+(27rn/Tp 2 )sin(2 1 rnT D /Tp)) 

+ (l/T p T 1 » 



2sin(2irnT D /Tp)/ (2irn) 


(A 4-c) 


X o = (E 1 -E 4 )T D /T p -(R 1 +R 2 )K 1 T D /Tp-((R 1 +R 2 )T 1 K 2 /Tp)(l-exp(-T D /T 1 )) 


(A 4-d) 


’ EaCTp-V/Tp 



2 (1-cos (2irnT D /T p ) ) / (2irn) 


(A 4-e) 


* rj 

Y 2 = (2/Tp) J exp(-t/T^)sin(nu> o t)dt * 

0 

(2/((l/T 1 ) 2 +(n«, 6 ) 2 ))([exp(-T I) /T 1 ))((-VCrpT 1 ))sin(2 1 rnT D /Tp) 
- (2trn/T p 2 )cos (2irnT D /T p ) ) + 2irn/T p 2 ) 


i 


(A 4-f) 


-43- 


T 

Y, - <-2/T p ) f P sin(n» 0 t)dt - 2(l-cos(2.nT D /T p ))/ (2*n) - Y p . (A 4-g) 
t D 

Let 

Z 1 * (E 1 -E 4 )(X 1 )-(R 1 +R 2 )(K 1 X 1 +K 2 X 2 )+E 3 X 1 (A 5-a) 


Z 2 = (E 1 -E 4 )(Y 1 )-(R 1 +R 2 )(K 1 Y 1 +K 2 Y 2 )+E 3 Y 1 . (A 5-b) 

The amplitude and phase angle of the n C ^ harmonic output 
voltage can be expressed in terms of and Z 2 as follows: 


V n = (( z i> 2+ ( z 2> 2 ) % and / v n " -tan" 1 (Z 2 /Z 1 ). (A 6) 


Hence the Fourier series for the output voltage is 


oo 

e 2 (t) = X Q + z^ V nQ cos(nu) o t-tan“ 1 (Z 2 /Z 1 )) . 


(A 7) 


Output Current i 2 



+ 


r Tp 

(2/T p ) J (K 3 exp(-(t-T D )/T 2 ))cos(nu, 0 
T D 


t)dt 


(A 8-b) 
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Again, it is convenient to evaluate the following set of 
integrals because they occur in several places in Equation (A 8) . 

T 

Y 4 = ((2exp(T D /T 2 ))/Tp) J P (exp(-t/T 2 >)cos(naj 0 t)dt - 

T D 

((2exp(T D /T 2 ))/((l/T 2 ) 2 +(na. o ) 2 ))((-I/Cr p T 2 ))exp(-Tp/T 2 ) 

- (exp ( -T D /T 2 ) ) ( (-1/ (T p T 2 ) ) cos (27rnT D /T p ) 

+ (27rn/Tp 2 )sin(2 1 mT D /T p )) < A 9-a) 


= ((2exp(Tjj/T 2 ))/Tp)/^exp(-t/T 2 )sin(nai t)dt ■ 

T D 

( (2 exp (T D /T 2 ) ) / ( (1/T 2 ) 2 +(na. 0 ) 2 ) ) ( (-2 1 m/T p 2 ) exp (-T p /T 2 ) 

- (exp (-T d /T 2 ) ) ( (-1/ (T p T 2 ) )sin(27rnT D /T p ) 

- (2Tm/Tp 2 )cos(27rnT D /Tp))) 


(A 9-b) 
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Y Q -K 1 T D /Tp+(K 2 T 1 /Tp)( 1 -exp(-T D /T 1 ))+(K 3 T 2 /Tp)(l-exp(-(Tp-T D )/T 2 )). (A 9-c) 


Let 

Z 3 - + K 2 X 2 + K 3 X 4 (A 10 -a) 

Z 4 " K 1 Y 1 + K 2 Y 2 + K 3 Y 5- (A 10_b) 

The amplitude and phase angle of the nth harmonic output 
current can be expressed in terms of Z 3 and Z^. 

! n - ((Z 3 ) 2 + <Z 4 ) 2 > % and/l n - -tan” 1 (Z 4 /Z 3 ) (A 11) 


Hence the Fourier series for the output current is 

i 2 (t) - Y + z I cos(nu) t - tan _1 (Z, /Z 3 )) . (A 12) 

n-1 o H 

At the input to the simple DC/DC chopper system, the 
instantaneous voltage and current is 

(A 13-a) 
(A 13-b) 
(A 13-c) 
(A 13-d) 

The in-phase, quadrature, and average Fouries series 
components for the input voltage are 




e l - E i - (K 1 +K 2 exp (-t/T L )); 0 < t <_ T D 


e l " E 1 


» T D - t - T P 


i l * Ki+K2 ex P(" t / T l) 


s °i t i T D 


i^ * 0 


i t d i c i T P 
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Input Voltage e^ 


T 

a n - (2/Tp) | D (E 1 -R 1 (K 1 +K 2 exp(-t/T 1 )))cos(nu, 0 t)dt 


+ (2E]/T p ) f cos(na) 0 t)dt 


(A 14-a) 


b n - (2/Tp) J (E 1 -R 1 (K 1 +K 2 exp(-t/T 1 )))sin(nu> 0 t)dt 
0 

T 

+ (2E 1 /T p )/ P sin(nai 0 t)dt 
T D 


(A 14 -b) 


T D T P 

(1/Tp) J (E i -R 1 (K 1 +K 2 exp(-t/T 1 )))dt + (E^Tp) J* dt . (A 14-c) 


Since the mathematical form of Equation (A 14) is similar 
to Equation (A 3), the amplitude and phase of the n c ^ 
harmonic is given by 


W x - -Rj.CKjXj. + K 2 X 2 ) 


(A 15-a) 


W 2 - -RxCK^ + K 2 Y 2 ) 


(A 15-b) 


W o aE l T D /T P'¥l T D /T p- ( R i T l K 2 /T p) (l-exp(-T D /T 1 ))+E 1 (Tp-T D )/T p 

(A 15-c) 


v n - ((W 1 ) 2 + (W 2 ) 2 ) % and/v n - -tan -1 (W^) 

i i 


(A 16) 


The Fourier series for the input voltage e^ is 

oo 

e., - + i cos(naj^t-tan -1 (W 2 /W^)) 


n-1 n i 


(A 17) 
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Input Current 1 ^ 

The in-phase, quadrature, and average Fourier series 
components for the input voltage are 

T 

a n “ < 2/T p> J' D (Ki+K 2 exp(-t/T 1 ))cos(n Uo t)dt (A 18-a) 

0 


r D 

b n «(2/T p ) J u (K 1 +K 2 exp(-t/T 1 ))sin(n«, 0 t)dt 


(A 18-b) 


a Q - a /T p) J (K 1 +K 2 exp(-t/T 1 ))dt. 

o’" 


(A 18-c) 


Since the mathematical form of Equation (A 18) is similar 
to Equation (A 8) , the amplitude and phase of the harmonic 
are given by 


W 3 - KjXj^ + K 2 X 2 


W 4 - K 1 Y 1 + K 2 Y 2 


W oo “ K 1 T D /T P + ^Ti/TpHi-expC-VT,)) 


D' *1' 


(A 19-a) 
(A 19-b) 
(A 19-c) 


I n i " ((W 3> 2 + < W 4> 2 ) % and An. “ -tan _1 (W 4 /W 3 ), 


(A 20) 


The Fourier series for the input current i^ is 


- 1 , 


i ■ W QO + Z I cos(na) Q t-tan (W 4 /W 3 )) 
n»l i 


(A 21) 
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In Section III of this report, Equations (35) and (36) 
represent the total average power input and output of a 
simple DC/ DC chopper system. The total average input and 
output power can be determined from Equations (24), (25), 

(26), (27), and (28) by multiplying e^ by i^ and e 2 by i 2 , 
respectively; and averaging the results over one period. 

The instantaneous input and output power is 

Pi - (E 1 -R 1 i 1 )i 1 ; 0<t<T D (A 22-a) 

p^ ■ 0 ; Tp <_ t <_ Tp (A 22-b) 


p o " (E 1 -E 4 -(R 1 +R 2 )i 1 )i 1 ; 0 < c < Tj 




3 2 


* T n i c 1 


(A 22-c) 
(A 22-d) 


The average total input and output power is 

T T 

P i * (E l /r p ) f Di l dt - (Ri/Tp) f^^dt 
J 0 0 


(A 23-a) 


T T Tp 

P 0 -((E 1 -E 4 )/T p ) J* D i 1 dt-((R 1 +R 2 )/T p ) j^i^dt- (E 3 /T p ) jT ^dt,^ 23 " b) 


respectively. Evaluating Equation (23) the results are 

P i " E l I o 1 " R 1 k 1 2 V T P ’ (2K 1 K 2 T 1 R 1 /T p )(l-exp(-T D /T 1 )) 

- (K 2 2 T 1 R 1 /(2T p ))(l-exp(-2T D /T 1 )) (A 24 -a) 


p 0 " ( E l- E 4> I o 1 ' (R 1 +R 2 )(K 1 2 T D /T p + (2K 1 K 2 T 1 /T p )(l-exp(-T D /T 1 ) 
+ (K 2 2 T ] /(2Tp))(l-exp(-2R D /T 1 )))-(E3K 3 T 2 /Tp) (1-exp (- (T p -T D )/T 2 ) ) 


(A 24 -b) 
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Where I is given by 

I 0 * (1/Tp) j ijdt - K^/Tp + (KjT^TpXl-expC-Tjj/T^) (A 25) 

1 o 

Equations (A 6), (A 7), (A 11), (A 12), (A 16), (A 17), 

(A 20), (A 21), (A 22), (A 23), (A 24), and (A 25) were pro- 
grammed on a computer. The program was arranged so that 
either the behavioral harmonic distribution up to and including 
the first 30 harmonics of the signals associated with simplified 
DC/DC chopper system could be studied as a function of on-time 
Tq and back emf, E 2 , or the power distribution over the same 
range of harmonics could be studied as a function of on-time 
and back emf. 
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DEFINITION OF SYMBOLS 

number of observed periods 

number of observed periods 

scanning parameter 

open- circuit voltage of the DC source (volts) 

shunt motor back emf (volts) 

semiconductor offset voltage (volts) 

semiconductor offset voltage (volts) 

instantaneous chopper input voltage (volts) 

instantaneous chopper output voltage (volts) 

frequency (hertz) 
corner frequency (hertz) 

fundamental frequency (hertz) 

transfer function in the current branch of the wattmeter 

transfer function in the voltage branch of the wattmeter 

amplitude of a sinusoidal current (amperes) 

amplitude of the n harmonic current (amperes) 

amplitude of the n harmonic chopper input current (amperes) 

phase angle of the n^ harmonic chopper input current (degrees) 

amplitude of the n c ** harmonic chopper output current (amperes) 

phase angle of the n^ harmonic chopper output current (degrees) 

DC component of current (amperes) 

instantaneous chopper input current (amperes) 

instantaneous chopper output current (amperes) 

minimum chopper output current (amperes) 

ratio of the 3-db frequencies 

deviation of the ratio of the 3-db frequencies 

shunt motor armature inductance (henrys) 
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N 

N 1 
NCF 
NPE 
NS PE 


n 

P APP 

^PP 

P 

°APP 

P(t) 

p(nT s ) 



data memory size 

highest harmonic of Interest 

normalized power correction factor 
formalized power error 
normalized sampled power error 
Integer 

approximate total average power (watts) 
approximate total average chopper Input power (watts) 

approximate total average chopper output power (watts) 

Instantaneous power (watts) 
sampled Instantaneous power (watts) 

equivalent resistance of the DC source and input leads (ohms) 

on-resistance of switch SW (ohms) . 

shunt motor armature resistance (ohms) 

chopper switch 

on-time of chopper (seconds) 

observation time (seconds) 

fundamental period (seconds) 

sampling time interval (seconds) 

electrical time constant of chopper system (seconds) 

electrical time constant of shunt motor (seconds) 

amplitude of a sinusoidal voltage (volts) 

amplitude of the harmonic voltage (volts) 

amplitude of the n c ^ harmonic chopper input voltage (volts) 

phase angle of the n C ^ harmonic chopper input voltage (degrees) 


L 



amplitude of the n c ** harmonic chopper output voltage (volts) 
phase angle of the n ttl harmonic chopper input voltage (degrees) 


DC component of voltage (volts) 
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x - normalized frequency (f/f^) 

9 ■ phase angle between sinusoidal voltage and current (degrees) 

* n - phase angle of the n th harmonic currant (degrees) 

* n ■ phase angle of the n c ^ harmonic voltage (degrees) 
w - angular frequency (radians/ sec. ) 

- comer angular frequency (radians/ sec.) 
u Q • fundamental angular frequency (radians/ sec. ) 
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FIGURE 1. General block diagram of a wattmeter. 
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figure 7. Maximum value of Che normalized sampled 

power error versus periods of observation. 
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figure 8. The mapiitude of che minimum value of 

the normalized sampled power error versus 
periods of observation. 





Ej - 100 volt 

Ej - 1 volt 

E^ - 1 volt 

Rj “ 0.05 ohm 

“ 0.03 ohm 

Rj ” 0.02 ohm 
Lj *0.0002 henny 


FIGURE 9. A schematic for a simple DC/DC chopper 



i 



-64 - 


•l 





— Fundamental component of the periodic signal. 

FIGURE 11. Typical operating characteristics for a DC /DC 
chopper with the DC component suppressed. 
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FIGURE 12. 


Typical input/output harmonic distribution for a 
simple DC/DC chopper (T n -0.0001 second and 

f? e 1 4 - - \ y 


E2“25 volts). 
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FIGURE 13-b. Input power measurement error versus number of 
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FIGURE I3-b. Input power measurement error versus number of 
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FIGURE 13-d. Input power measurement error versus number of 
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FIGURE 14-a. Flow chare for a microcomputer 
wattmeter. 
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FIGURE 14-b. Flow chare for a microcomputer 
wattmeter (continued) . 



FIGURE 14-c. Flow chare for a microcomputer 
wattmeter (continued). 
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figure 14 -e. Flow chart for a microcomputer 
wattmeter (continued) . 
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FIGURE 15-a. Microcomputer output, microcomputer NSPE, and 
theoretical NSPE veraua periods observed. 
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FIGURE 15-b. Microcomputer oucput, microcomputer NSPE, and 
theoretical NSPE versus periods observed. 
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FIGURE 15 -c. Microcooputtr output, nlerocoaputar NSPE, and 
tbaoratleal NSPE varsus pariods obstrvsd. 
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FIGURE 13-d. Microcomputer output, nlcroeoaputer NSPE, and 
thaoratlcal NSPE versus periods observed. 
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FIGURE 15-e. Mlcroconpucer output, nicrocoaputer NSPE, and 
theoretical NSPE veraua period* obaerved. 
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FIGURE 16. Mleroeonputer output versus the periods observed 
for different phase angles. 
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